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' Abstract. Some selected concepts for the solar activity cycle are briefly reviewed. Cycle mod- 
I ulations through a stochastic a effect are being identified with limited scale separation ratios. 
Three-dimensional turbulence simulations with helicity and shear are compared at two different 
scale separation ratios. In both cases the level of fluctuations shows relatively little variation 
with the dynamo cycle. Prospects for a shallow origin of sunspots are discussed in terms of the 
negative effective magnetic pressure instability. Tilt angles of bipolar active regions are discussed 
as a consequence of shear rather than the Coriolis force. 
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p/j 1. Solar cycle 

• The solar cycle manifests itself through spots at the Sun's surface. To understand 

' activity variations, we have to understand not only their source, but also the detailed 

O I connection between variations in the strength of the dynamo and its effect on the number 

Q and size of sunspots. In this paper, we address both aspects. 

^ ■ The physics of the solar cycle is not entirely clear. The models that work best are 

C/3 I not necessarily those that would emerge from first principles. Even the reason for the 
equatorward migration of the activity belts is not completely clear. Following Parker's 

' seminal paper of 1955, this migration seemed to be a simple property of an all dynamo, 

T-H ' i.e., a dynamo that works with a effect and shear. What matters for equatorward migra- 

^ tion is not the gradient in the latitudinal direction, but that in the radial one, dfl/dr. 
However, in the bulk of the convection zone, dVl/dr is mostly positive. This, together 

' with an a effect of po sitive sign in the northern hemisphere results in poleward migra- 



tion (jYoshimuralll975l ). which is not what is observed. On the other hand, according to 
the flux transport dynamos, magnetic fields are advected by the meridional circulation. 
Assuming that there is a coherent circulation with equatorward migration at the bottom 
of the convection zone, this would then turn the dynamo wave around so as to explain 
the solar butterfly diagram and that sunspots emerge from progressively lower latitudes 



( Choudhuri et al. 1995: Dikpati fc Charbonneau.l999j : iGuerrero fc de Gouveia Dal Pino 



2008() . This requires that most of the field resides at the bottom of the convection zone. 



' Moreover, the a effect is taken to be non-vanishing only near the very top of the convec- 
^ . tion zone, i.e., the mean electromotive force has to be written formally as a convolution 
" " ' of t he mean magnetic field with a n integral kernel to account for this non-locality (see. 



e.g.. [Brandenburg fc Kapvlall2007l ). Furthermore, from the observed tilt angles of bipolar 



regions it is inferred that the Sun's mag netic field at the bottom of the convection zone 



reaches strengths of the order of 100 kG (jD'Silva fc Choudhurilll993l ). which is nearly 100 



times over the equipartition value. Finally, there are assumptions about the turbulent 
magnetic diffusivity. In all cases, the magnetic diffusivity in the evolution equation for 
the toroidal field in the bulk of the convection zone is rather small, below 10^^ cm^ s^^ 
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see, e.g.. IChatteriee fc Choudhuril 120061 ) . The magnetic diffusivity for the poloidal field 



is assumed to be larger and similar to the values expected from mixing length theory 
(see below). 

In any case, these assumptions are hardly in agreement with standard formulae that 
the magnetic diffusivity is given by iru^^g, where t is the turnover time and Wrms is 
the rms value of the turbulent velocity. The turnover time is r = (wrms^f)"^: where fcf 
is the waven umber of the en ergy-carrying eddies. This result for rjt is well confirmed by 
simulations ( Sur et a02008l) . For the Sun, mixing length theory appears to be reasonably 
good and gives rjt ~ (1...3) x 10^^ cm^ s~^. Also, in contrast to the assumptions of some 
fiirx transport dynam o models, a strong degre e of anisotropy of the rj tensor is not 
expected from theory ()Brandenburg et al. 2011 ). 

An alternate approach is to use turbulent transport coefficients from theory, which 
give rise to what is called a distributed dynamo, i.e., the induction effects are non- 
vanishing and distributed over the entire convection zone. In addition, there is the hy- 
poth esis that the near-su rface shear layer may be important for the equatorward migra- 
tion ( Brandenburg! 2005 ). but this has never been confirmed by simulations either. In 
any case, based on such models one would not expect there to be a 100 kG magnetic 
field, but only a much weaker field of around 0.3-1 kG. This calls then for an alterna- 
tive explanation for the magnetic field concentrations of up to 3kG seen i n sun spots 
and active regions. Various proposals were already discussed in BrandenbureJ ( 20051 ). and 
meanwhile there are direct numerical simulations (DNS) confirming the validity of the 
physics assumed in one of those proposals. This will be addressed in Section [31 



2. Cycle modulation 

Early ideas for cycle modulations go back to Tavakol ( 19781 ) who argued that the solar 
cycle may be a chaotic att r actor. This expla nation became very popular in the follow- 
ing years ( RuzmaikinI 1 1 98 1 : Weiss et al. 1984). These ideas were elaborated upon in the 
framework of low-order truncations of mean-field dynamo models, having in mind that 
the same idea applies also to the underlying fully nonlinear three-dimensional equations 
of magnetohydrodynamics. Another line of thinking is that in mean-field theory (MFT) 
the physics of the cycle m odels can be expl ained b y ra ndom fluctu ations in the turbu- 
lent transport coefficients (Choudhuri et al. 1992t IMoss et al. I ll992t ISchmitt et aLlll996t 
Brandenburg fc Spiegell 2008). For high-dimensional attractors there is hardly any dif- 
ference between both approaches. A completely differe nt proposal for cy cle modulation 
is related to variations in the meridional circulation (jNandv et a/.ll201ll ). This proposal 
still lacks verification from DNS of a dynamo whose cycle period is indeed controlled by 
meridional circulation. By contrast, fluctuation s in the turbulent transp ort coefficients 
have indeed been borne out by simulations (see iBrandenburg et al. 20081 ) . 

To illustrate this, let us now consider a physical realization of a simple aVl dynamo in 
a periodic domain. In the language of MFT, this corresponds to solving the following set 
of mean-field equations. 



— = V X {U xB + £-rjnoJ) 



(2.1) 



in a Cartesian domain, (x, y, z), in one dimension, — tt < z < tt, where U — Us = 
(0,5^,0) is a linear shear flow velocity, J = V x i3///o is the mean current density, /io 
is the vacuum permeability, i] is the microphysical (molecular) magnetic diffusivity, and 



(2.2) 
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is the mean electromotive force. In DNS, on the other hand, one often solves the equations 
for an isothermal gas with constant sound speed Cg, 

f)7-$ 

— = Vx{UxB- rj^ioJ) , (2.3) 

together with corresponding equations governing the evolution of the turbulent velocity 
U. 

In the following we present results of simulations using shearing-periodic boundary 
conditions. To maintain the solenoidality of the magnetic field, we write B = 'V x A 
and solve for the magnetic vector potential A. Using in the following the velocity for the 
deviations from the shear flow, [/, our equations are 

8 A 

— + Us ■VA = -SAyX + U X B + T]V^-A, (2.4) 

^ = -SU^y -clV\np + f + -{JxB + V- 2vpS) , (2.5) 
Dt p 

where D/Dt = d/dt+ {U + Us) ■ V is the advective derivative with respect to the total 
flow, U + Us, p is the gas density, v is the viscosity, Sij — ^{Ui^j + Uj^i) — ^Sij^ -U is the 
trace-less rate of strain matrix, and / is a forcing function that drives both turbulence and 
a linear shear flow. Alternatively, turbulence can also be the result of some instability 
(Rayleigh-Benard instability, magneto- rotational instability, etc). In the following we 
restrict ourselves to a random forcing function with wavevectors whose modulus is in a 
narrow interval around an average wavenumber fcf. This has the additional advantage 
that we can arrange the forcing function such that it has a part that is fully helical, i.e., 
V X / = fcf/ (the part driving the shear flow is of course non- helical) . Because of the 
presence of helicity, we should expect there to be an a effect operating in the system, but 
if the number of turbulent eddies in the domain is not very large, there can be significant 
fluctuations in the resulting a effect. 

Important control parameters are the magnetic Reynolds and Prandtl numbers, Kcm = 
Wrms/jyfcf and PrM — v/rj. In addition, there is the shear parameter defined here as 
Sh = The smallest possible wavenumber in a triply-periodic domain of size 

L X L X L is ki — 2'k/L. For the purpose of presenting exploratory results, we restrict 
ourselves here to a resolution of 64'^ meshpoints. We use the fully compressible Pencil 
CodeQ for all our calculations. 

In Figure [T] we present the results of two simulations with scale separation ratios fcf/fci 
of 1.5 and 2.2. In both cases, fcf/fci is still relatively small, but the difference in the 
results is already quite dramatic. For fcf/fci = 2.2 the cycle is more regular while for 1.5 
it is quite erratic. We show the toroidal field (i.e. the component By in the direction of 
the mean shear fiow) at an arbitrarily chosen mesh point as well as its squared value 
(which could be taken as a proxy of the sunspot number), the mean magnetic energy 
in the full domain, as well as its contributions from the mean and fiuctuating fields. In 
all cases, the magnetic field is normalized by the equipartition value, Bcq = ^//ioPoUrms , 
where po = (p) is the mean density, which is conserved for periodic and shearing-periodic 
boundary conditions. Here, angular brackets denote volume averages and overbars are 



t http://www.pencil-code.googlecode.com 
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Figure 1. Time sequences of By, By, {B^) (solid line) together with (B^) (dotted lino), and 
(b^). The field is always normalized by Beq- Here, Rcm = 22 for fef/fci = 1.5 (left column) 
and Rcm = 9 for fcf/fci = 2.2 (right column). In both cases, PrM = 5. The shear parameter is 
Sh ~ —2 in both cases. 

defined as xy averages, so 

B{z, t) = j B{x, y, z, t) dx dy/L^Ly, (2.7) 

which implies that 

{B^) = {W) + {b''). (2.8) 
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Figure 2. Time series of By/Bcq for fcf/fci — 1.5 (left column) and fcf/fci — 2.2 (right column) 
for different values of Sh. Again, ReAf — 22 for fcf/fci = 1.5 (left column) and Ksm = 9 for 
fcf/fci = 2.2 (right column), and PrM = 5 in both cases. For fcf/fci = 2.2 the oscillations tend to 
become less sinusoidal for larger values of Sh, while for fcf/fci = 1.5 there are strong fluctuations 
that tend to become somewhat weaker for larger values of Sh. 

One sees that (b^) shows fluctuations that are not strongly correlated with the varia- 
tions of the mean field (Figure [T]) . This is important because the lack of a correlation is 
sometimes used to argue that the Sun's smah-scale magnetic field must be created by a 
local small-scale dynamo and disconnected with the large-scale dynamo. 

There are several other interesting differences between the two cases. The cycle period 
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is given by Wcyc ~ Vtkf i Kapvla fc Brandenburd2009l ). and with rji « Urms/3fcf ( Sur et al 



.. „ '^cyc ^ , , , ,, , . , 

12008,) we have Ucyc ~ jWnnsfcf (^i/^f)^; thus the normahzed cycle period is TcycWrms^f 
27rurmsfcf/wcyc ~ 67r(fcf/fci)^ « 91 for fct/fci = 2.2, which agrees with the result shown in 
the upper right panel of Figure [T] Next, for kf/ki — 2.2 the mean magnetic energy and 
{B'^) are about 4 times larger than for fcf/fci = 1.5. This value is larger than the one 
expected fro m the theory where this ratio sh ould be equal to the ratio of the respective 
values of k[ ( Blackman fc Brandenburgll2002h . namely 2.2/1.5 ~ 1.5. 



Next, we ask how the properties of the dynamo change as it becomes more supercritical. 
This is shown in Figure [5] where we plot time series of By/Beq for fcf/fci = 1.5 and 
k[/ki — 2.2 for different values of Sh. Note that for fcf/fci = 2.2 the oscillations become 
less sinusoidal for larger values of Sh, while for k{/ki = 1.5 there are strong fluctuations 
that become somewhat weaker for larger values of Sh. 

A more quantitative way of assessing the properties of the large-scale dynamo is by 
looking at the scaling of the magnetic energy of the mean field and the cycle frequency 
as a function of the nominal dynamo number, D — CaCs, where 

C„ = a/r/Tfci = teffcf/fci, Cs ^ S/rirki ^ 3iSh{kf/ki)^, (2.9) 

are non-dimensional numbers measuring the expected value of the a effect (assuming 
a « ^t{u; ■ It), rjt « |t(m^), with r = (wrmsfcf)"^ and (a; • m) « fct(it^)), as well as the 
shear or fi effect. In these approximations, rj^ = rj + rjt is the expected total magnetic 
diffusivity and t = (1 -I- 3/Rej\/)~^ is a correction factor that takes into account finite 
conductivity effects resulting from the fact that rjt ^ Vt- The results are shown in 
Figure [H 

For fcf/fci = 1.5, the scaling of {B'^)/Bl with D suggests that t he critical value is be- 



tween 1 and 2, i.e., somewhat smaller than the theoretical value of 2 ( Brandenburg fc Subramanian 



2005h . For fcf/fci = 2.2 the critical value is < 1. The cycle frequencies are approximately 



independent of D, except that for /cf/fci = 1.5 there is a sharp drop for D > 10. Owing to 
fluctuations, a Fourier spectrum of the time series is not sharp but has a certain width. 
We determine the quality or width, wq, by fitting the spectrum to a Gaussian propor- 
tional to P{ijj) ^ exp[— (w — uJof /2'Wq]. Also the values of wq, shown in the last panel 
of Figure [3l are approximately independent of D. For fcf/fci = 2.2, wq is substantially 
smaller than for fcf/fci = 1.5. This indicates that the cycle period is better defined for 
larger scale separation ratios. 



3. Active regions and their inclination angle 

One should expect that the sunspot number depends in a complicated way on the 
magnetic field strength. If sunspots are indeed relatively shallow phenomena, the field 
must be locally concentrated to field strengths of up to 3 kG. A candidate for a mechanism 
that can concentrate mean fields of ^ 300 G, which is about 10% of the local equipartition 
field strength, is the negative effective magnetic pressure instability (NEMPI). This is 
a remarkable phenomenon resulting from the suppression of turbulent pressure by a 
moderately strong large-scale magnetic field. This suppression is stronger than the added 
magnetic pressure from the mean field itself, so the net effect is a negative one. 

The fact that this phenom enon can lead to an instability in a stratified layer wa s first 
found in mean- field m odels (iBrandenburg eraD l2010l l2011bt iKapvla et a]]l201ll ). and 
more recently in DNS ([Brandenburg et aLll2011ai r However, NEMPI has not yet been 
able to explain flux concentration in the direction along the mean magnetic field, i.e., the 
large-scale structures remain essentially axisymmetric. To discuss the theoretical origin 
of this, we need to look at the underlying mean-field theory. Similar to the effective 
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Figure 3. Scaling of the relative magnetic energy of the mean field, {B'^)/B'^q, normalized 
cycle frequency, and normalized quality versus dynamo D = CaCs- 



magnetic diffusivity in the mean electromotive force, the sum of Reynolds and Maxwell 
stresses from the small-scale field depends on the mean magnetic field in a way that looks 
like a Maxwell stress from the mean field, but with renormalized coefficients. The concept 
of expressing the Reynolds stress from the fluctuating velocities, UiU], by the mean flow 
U is of course familiar and leads to the usual turbulent viscosity term, —i^t{Ui,j + Ujs)- 
However, in the presence of a mean magnetic field, symmetry arguments allow one to 
write down additional components, in particular those proportional to SijB"^ and BiBj. 
The sum of Reynolds and Maxwell stresses from the fluctuating velocity and magnetic 
fields is given by 

= plMU] - bibj/flo + ^b'^/fJ'O, (3.1) 
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where the superscript f indicates contributions from the flu ctuating fiel d. Expres sing If,,,- 
in terms of the mean field, t he leading terms are ( Kleeorin et al..l990i : iKleeorin fc: Rogachevski: 
1994 iKleeorin et alJIlQOa iRogachevskii fc Kleeorinll2007l) 



— f 



(3.2) 



where the dots indicate the presence of additional terms that enter when the effects 
of stratification affect the anisotropy of the turbulence further. Note in particular the 
definition of the signs of the terms involving the functions qs{B) and qp{B). This becomes 
obvious when writing down the mean Maxwell stress resulting from both mean and 
fluctuating fields, i.e.. 



B,Bj/^io 



= -(1 - q^)B,Bjlno + \{l - qp)S^jB'^/^io 



(3.3) 



A broad range of different DNS have now confirmed that qp is positive for Re^f > 1, but 
qs is small and negative. A positive v alue of qs (but with large error bars) was originally 
reported for unstratified turbulence (iBrandenburg et al. 2010 ). Later , stratified simula- 
tions with isothermal s table stratification ( Brandenburg et al. 2011bf ) and convectively 
unstable stratification (jKapvla et aZ.1 120111 ) show that it is small and negative. Never- 



theless, qp{B) is consistently positive provided Re a/ > 1 and B/Bcq is below a certain 
critical value that is around 0.5. This implies that it is probably not possible to pro- 
duce flux concentrations stronger than half the equipartition field strength. So, making 
sunspots with this mechanism alone is maybe unlikely. 

The significance of a positive gs value comes from mean-field simulations with gs > 
indicating the formation of three-dimensional (non-axisymmetric) flux concentrations 
( Brandenburg al. 2010l). This result was later identified to be a direct consequence of 
having gs > ( Kemel et al.ll2011 ). Before making any further conclusions, it is important 
to assess the effect of other terms that have been neglected. Two of them are related to 
the vertical stratification, i.e. terms proportional to QiQj and ,g,; i?, -|- q^Bj with q b eing 
gravity. The coefficient of the former term seems to be small ([Kapvla et g/.l 120111 ) and 
the second only has an effect when there is a vertical imposed field. However, there could 
be other terms such as JiJ j as well as JiBj and J jBi that have not yet been looked at. 

Yet another alternative for causing flux concentrations is the suppression of turbulent 
(convective) heat transport whi ch might even be strong enough to explain the formation 
of sunspots (jKitchatinov fc Mazur. 2000.) . It is conceivable that effects form heat transport 
become more important near the surface, so that a combination of the negative effective 
magnetic pressure and the suppression of turbulent heat transport are needed. Another 
advantage of the latter is that this mechanism works for vertical fields and is isotropic with 
respect to the horizontal plane, so one should expect the formation of three-dimensional 
non-axisymmetric structures. 

Once an active region forms and it is bipolar, we must ask ourselves how to explain the 
observed tilt angle. This question cannot be answered within the framework of NEMPI 
alone, but i t requires a c onnection with the underlying dynamo. Here, we can refer to 
the work of lBrandenbrirel [2005) where bipolar regions occur occasionally at the surface 
of a domain in which shear-driven turbulent dynamo action was found to operate; see 
Figure m The reason for the tilt is he r e not the Coriolis force, as is usujilly assum ed, bu t 
shear; see also iKosovichev fc Stenflol ( 20081 ). In the simulations of iBrandenburgl ( 2005 ). 
shear was admittedly rather strong compared with the turbulent velocity, so the effect is 
exaggerated compared to what we should expect to happen in the Sun. However, even 
then there are a few other problems. One of them is that the bipolar regions appear 
usually quite far away from each other (Figure |4|). This may not be realistic. On the 






Figure 4. Magnetograms of the radial field at the outer surface on the northern hemisphere 
at different times for a simulation presented in iBrandenburd (|2005l 'l. Light shades correspond to 
field vectors pointing out of the domain, and dark shades correspond to vectors pointing into 
the domain. The elongated rings highlight the positions of bipolar regions. Note the clockwise 
tilt relative to the y (or toroidal) direction, and the systematic sequence of polarities (white left 
and dark right) corresponding to By > 0. Here, the z direction corresponds to latitude. 



other hand, it is not clear how to scale this model to the Sun. In this model, the scale 
separation ratio is rather small, so the extent of the bipolar regions is comparable to a 
few times the convective eddy size, which, for the solar surface, is not very big (a few 
Mm). However, on that small scale one would not expect the effects of differential shear 
to be very important. 

Another issue is that in the model of Brandenburgl ( 2005 ). bipolar regions occur only 
occasionally. To illustrate this, we show in Figure[S]the resulting magnetograms for three 
times that are separated by about 2 turnover times. Clearly, other structures can appear 
too and the field is not always bipolar. 
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Figure 5. Similar to Figure [4^, showing also the results 4 and 2 turnover times earlier. The 

last panel is identical to Figure |4^. 



4. Conclusions 

It is clear that the magnetohydrodynamics of mean magnetic and velocity fields is 
quite rich and full of important effects. The standard idea t hat sunspots and bipolar 
regions form as a result of an instabili ty in the tachocline ( Gilman fc Dikpatil 20001 : 
Gaily et al. l l2003l:IParfrev fc Menoul[2007h may need to be re-examined in view of several 
new alternative proposals being on the horizon. In addition to comparing model with 
observations at the solar surface, there are ways of comparison b oth beneath the surf ace 
and above. Particularly exciting are the recent determinations of lllonidis et ali (j201l[ ) of 
some sort of activity at w 60 Mm depth. It is also of interest to explain magnetic activity 
in the solar wind, and especially its magnetic helicity which has rec ently been found to 



be bi -helical, i.e., of opposite signs at large and small length scales ([Brandenburg et al 



201 If ) and positive at small length scales in the north. This is particularly interesting, 
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be cause such a result ha s recently been reproduced by distributed dynamo simulations 



of IWarnecke et al\ (|201lh who also find positive magnetic helicity at small length scales 
in the north. More detailed and varied comparisons between the different approaches are 
thus required because we can more fully understand the Sun's activity cycles and their 
long term variations. 
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Discussion 

Munoz-Jaramillo: If sunspots are produced at the surface, which processes would 
lead to their decay? 

Brandenburg: I think it could be the continued submersion of magnetic structures. 
The system remains time-dependent and new structures will form near the surface, while 
old ones disappear from view. 

Luhmann: Axel, you are one of the few dynamo modelers that include the corona and 
larger heliosphere in your models and thinking. How important is that to the results of 
the models, and do models not including that aspect have compromised results? 

Brandenburg: The magnetic helicity flux divergence is crucial for alleviating catas- 
trophic alpha quenching; see the next talk by Candaleresi et al. (2011). Helicity fluxes 
through surface carry about 30of the total; the rest goes through the equator. The ob- 
served magnetic helicity spectra support our understanding in terms of the magnetic 
helicity evolution equation. Models not including helicity fluxes suffer artificially strong 
catastrophic quenching, but only if their magnetic Reynolds numbers are really large. 

Choudhuri: Flux rise simulations based on the idea that the toroidal field forms in the 
tachocline explained Joy's law other characteristics of sunspot graphs. Can these results 
be recovered if sunspots form from near-surface fields? 

Brandenburg: The tilt angles of near-surface produced flux concentrations is deter- 
mined by latitudinal shear, as was demonstrated by Brandenburg (2005). 

Nandi: For the near-surface, shallow dynamo to work, you need processes that are slow 
enough to store, amplify, and transport fields on 11 year timescales. However, in the 
upper convection zone, eddy turnover timescales are short. Can you comment on how 
you reconcile this with your shallow dynamo? 

Brandenburg: Magnetic helicity conservation is generally responsible for prolonging 

the time scales. In fact, the partial alleviation of catastrophic quenching by magnetic heli- 
city fluxes means that the timescales are not infinite. The dynamo period is proportional 
to {adfl/dr)~^/^ , so a quenching prolongs the period. 



